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ABSTRACT 


This  report  contains  simple  formulas  and  graphs  that  can  be  used 
to  obtain  quick  estimates  of  the  total  current  induced  in  a  buried 
cable  by  an  incident  plane-wave  exponential  pulse.  Formulas  and  graphs 
are  also  given  for  estimating  the  internal  voltage  and  current  in  cables 
with  tubular  shields  when  the  total  current  in  the  cable  is  an  exponen¬ 
tial  pulse,  and  when  the  current  is  that  induced  in  the  cable  by  the 
incident  exponential  pulse. 
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This  report  has  been  prepared  to  assist  the  worker  in  EMP  coupling 
studies  in  obtaining  quick  estimates  of  the  current  induced  in  buried 
cables  and  of  the  internal  current  and  voltage  induced  in  shielded 
cables.  An  attempt  has  been  made  to  organize  the  material  in  a  format 
that  permits  the  desired  result  to  be  obtained  with  a  minimum  of  addi¬ 
tional  calculation  and  without  extraneous  development  of  the  theory 
involved.  Following  the  presentation  of  the  results,  however,  a  brief 
explanation  of  their  source  is  provided  to  give  the  interested  reader 
more  insight  into  the  nature  and  quality  of  the  results.  It  is  felt 
that  such  information  will  be  useful  to  the  user  in  establishing  the 
applicability  and  limitations  of  the  results,  although  the  range  of 
applicability  is  summarized  at  the  outset  in  each  case  under  the  heading 
"Class  of  Problems." 

The  bulk  (total)  current  in  the  buried  cable  is  presented  for  a 
plane-wave  incident  on  the  surface  of  the  earth.  The  incident  waveform 
is  an  exponential  pulse;  however,  the  decay  time  constant  is  arbitrary 
(within  limits) ,  so  that  incident  field  waveforms  approaching  an  impulse 
(very  short  decay)  or  a  step  function  (very  long  decay)  can  be  obtained. 
The  other  parameters  that  affect  the  cable  response  are  the  direction 
of  arrival  and  the  soil  conductivity.  The  direction  of  arrival  is  a 
geometric  function  that  is  a  separable  directivity  coefficient  D(\[f,cp) 
of  the  current  spectrum  or  waveform.  The  soil  conductivity  is  incorpo¬ 
rated  through  a  soil  time  constant,  t  =  e  /a ,  for  convenient  normaliza- 

e  o 

tion  of  the  cable  responses.  Thus  the  bulk  current  in  the  cable  is 
defined  in  terms  of  the  three  parameters,  t  (the  incident-pulse  decay 
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time  constant),  t  (the  soil  time  constant),  and  D(i)f,cp)  (the  directivity 
e 

function) ,  in  addition  to  the  amplitude  of  the  incident  field.  A  single 
waveform  in  normalized  time  t/T  describes  the  bulk  current  in  a  long 
buried  cable. 

The  bulk  current  results  are  given  only  for  the  case  in  which  the 
ground  behaves  as  a  conductor.  For  most  angles  of  incidence  and  soil 
conductivities,  this  limitation  is  inconsequential  for  the  EMP  spectrum 
because  the  relative  magnitude  of  the  high-frequency  spectrum  where  the 
soil  behaves  as  a  lossy  dielectric  is  small  and  the  attenuation  of  the 
field  for  typical  cable  burial  depths  at  these  frequencies  is  significant. 
Furthermore,  the  cable  shield  usually  attenuates  these  frequencies  so 
severely  that  they  are  completely  negligible  in  the  internal-voltage  and 
-current  analysis. 

The  internal  voltage  and  current  are  presented  for  cables  with 

uniform  tubular  shields.  The  important  cable  parameters  affecting  these 

internal  responses  are  the  dc  resistance  per  unit  length  of  the  shield, 

2 

the  length  of  the  cable,  the  shield  diffusion  time  constant  Tg  =  juaT  , 
and,  for  the  core  current,  the  characteristic  impedance  of  the  core-to- 
shield  transmission  line  ("core"  is  used  to  designate  all  the  conductors 
covered  by  the  shield) .  The  internal  responses  are  presented  for  an 
exponential  pulse  of  current  in  the  shield.  In  addition,  for  electrically 
short  cables  the  internal  responses  of  a  buried  cable  to  an  incident 
plane-wave  exponential  pulse  are  presented.  Only  open-circuit  or  matched 
(characteristic  impedance)  terminations  are  used  for  the  internal 
responses  presented  here. 

Finally,  in  an  effort  to  make  this  document  a  more  useful  reference 
for  the  cable-coupling  analyst,  appendices  containing  transmission-line 
relations,  electromagnetic  properties  of  common  shielding  materials,  a 
wire  table,  and  a  soil-conductivity  map  have  been  included.  Appendix  A 
is  a  list  of  the  symbols  used,  and  their  meanings. 
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II  BULK  CURRENTS  IN  BURIED  CABLES 


A.  General 

This  section  presents  formulas  for  calculating  the  total  current 
induced  in  a  buried  cable  by  an  incident  uniform  plane  wave,  for  bare 
and  thinly  insulated  cables.  The  incident  electromagnetic  pulse  is 
assumed  to  have  an  exponential  waveform  Eq  exp  (-t/T).  The  wave  is 
assumed  to  arrive  on  the  surface  of  the  earth  from  a  direction  defined 
by  an  elevation  angle  if  and  an  azimuth  angle  cp  as  illustrated  in 
Figure  1.  The  depth  of  burial  of  the  cable  is  small  compared  to  a  skin 
depth  in  the  soil,  so  that  the  fields  at  the  cable  depth  are  essentially 
the  same  as  those  at  the  surface. 

DIRECTION  OF  PROPAGATION 


FIGURE  1  ILLUSTRATION  OF  COORDINATES  FOR  CABLE  AND  ANGLE  OF  ARRIVAL 
OF  TRANSIENT  WAVE 
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FIGURE  2  SKIN  DEPTH  8  IN  SOIL  AS  A  FUNCTION  OF  FREQUENCY  AND  SOIL  CONDUCTIVITY 


In  all  the  bulk-current  formulas,  the  soil  is  considered  to  be  a 
good  conductor,  as  defined  by  a  »  cue.  The  wave-propagation  factor  in 
the  soil,  and  the  propagation  factor  for  bare  conductors,  is  thus 


-J- 


(a  + 
o 


ju)e) 


jcuM'  a 
o 


1±J 

6 


(1) 


where  6  is  the  skin  depth  in  the  soil.  A  plot  of  skin  depth  as  a 
function  of  frequency  and  soil  conductivity  is  given  in  Figure  2.  The 
dashed  line  along  the  right-hand  side  of  Figure  2  indicates  the  limit  of 
validity  of  the  approximation  y  »  (1  +  j)/6.  Although  this  limit  falls 
between  a  few  hundred  kHz  and  a  few  MHz  for  typical  soil  conductivities, 
most  cable  shields  have  very  large  attenuation  at  frequencies  above 
100  kHz,  so  that  the  bulk-current  spectrum  at  these  frequencies  is 
usually  of  secondary  interest. 

The  current  induced  in  a  buried  cable  is  computed  by  treating  the 

cable  as  a  transmission  line  that  is  driven  by  the  component  of  the 

underground  electric  field  that  is  parallel  to  the  axis  of  the  cable. 

As  illustrated  in  Figure  1,  the  z  axis  is  parallel  to  the  axis  of  the 

cable,  and  the  direction  of  arrival  of  the  incident  wave  will  be 

specified  by  the  elevation  angle  \|r  and  the  azimuth  angle  cp.  The  phase 

of  the  incident  wavefront  at  the  air/earth  interface  is  of  the  form 

k  cos  ijr  cos  cp,  where  k  =  (u/(i  €  ,  and  the  component  of  the  electric 

o  o 

field  parallel  to  the  cable  will  be 

E  =  E(cu)e  Z(1  -f  R  )  sin  cp  (2) 

z  h 


for  horizontal  polarization,  and 


p-  2 

E  =  E(tu)  e  (1  -  R  )  cos  cp  sin  \Ir 

z  V 


(3) 
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for  vertical  polarization,  where 

k/  =  k  cos  t|r  cos  cp 

E(cju)  =  Incident  field  strength 

R  =  Reflection  coefficient  for  the  air/earth  interface  when  the 
h  incident  wave  is  horizontally  polarized 

R  =  Reflection  factor  when  the  incident  wave  is  vertically 
v 

polarized. 


These  reflection  factors  are  given  by 


si 


Rh  = 


n  ft  -  J e  (1  +  -2—)  - 
y  r  joue 


2 

cos  l|f 


sin  ft  +  e  (1  +  - 

Y  r  jtue 


2 

cos  \|f 


(4) 


and 


R 


e  (1  +  - )  sin  Ur 

r  joue _ 


4 


a  2 

e  (1  +  - )  -  cos  Ur 

r  joue _ 


v  a  / a”  2^ 

e  (1  +  - )  sin  *  +  / e  (1  +  “ — )  -  cos  i|r 

r  juue  V  r  jcoe 


(5) 


where  a  is  the  soil  conductivity,  e  is  the  dielectric  constant  of  the 


soil,  and  e  =  e  e  .  For  cr  »oue,  the  terms  (1  ±  R)  reduce  to 
o  r 


1  +  R  « 


2  sin  ft 


h  J a 

v  ju)e0 


(6) 


1  -  R  <=« 


v  rm 
VJujec 


sin  t|r 


(7) 


The  current  induced  in  the  cable  is  then 


I  (z,uj)  =  [K  +  P(z)]e"YZ  +  [K  +  Q(z)  ] 
1  £ 


yz 


(8) 
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where  y  =  /ZY,  the  propagation  factor  for  the  cable,  with  the  soil  as 

the  return  conductor.  K  and  K  are  constants  given  by 

1  ^ 


Y  2  Y  2 
YZ1  P2P(Z2)e  "  Q(zi)e 


Y<VV 


-  p  p  e 
±k2 


y(z2-z1) 


K2  =  P26 


Yzi  -Yzi 

-yz  p  Q(z  )e  -  P(z  )e 
2X1  2 


y(VV 


P  P„e 
K1  2 


-7(^-2,) 


in  which  the  reflection  coefficients  p  and  p  are  given  by 

1  2 


z  -  z  z  -  z 

JL _ 0  _  _2 _ 0 

P1  -  Z  +  Z  ’  P2  ~  Z  +  Z  * 
10  2  0 


and  the  source  functions  P(z)  and  Q(z)  are 


P(z)  =  ~  J*  eYV  E  (v)  dv 


0  z. 


— -  J  e  Y  E  (v)  dv 
2Z  J  z 

0  z 


The  impedances  Z  and  Z  are  those  that  terminate  the  cable  shield 
1  2 

at  the  ends  z  and  z  ,  respectively.  When  Z  =  Z  =  Z  ,  the  values  of 
12  12  0 

p  ,  p  ,  K  ,  and  K  are  all  zero.  This  is  the  case  of  infinitely  long 
12  1  2 

cables  or  cable  shields  terminated  in  their  characteristic  impedances. 


7 


The  characteristic  impedance  is  obtained  from  the  impedance  per 
unit  length,  Z,  and  the  admittance  per  unit  length,  Y,  from  Zq  =  / Z/Y . 

B.  Current  Far  From  the  Ends  of  a  Long  Cable 

1 .  Class  of  Problems 

The  response  characteristics  presented  here  apply  to  points 
far  (several  soil  skin  depths)  from  the  ends  of  long  buried  cables  that 
have  no  insulation  (or  have  insulation  that  is  thin  compared  to  the  cable 
rsdius )  over  the  sheath  (shield  or  armor)*  These  response  character 
istics  also  apply  to  cables  with  semiconducting  jackets  so  that  the 
metal  sheath  is  effectively  in  contact  with  the  soil.  The  responses 
are  based  on  the  assumption  that  the  soil  behaves  as  a  conductor 
(a  >  cue),  and  that  the  radius  of  the  cable  is  small  compared  to  a  skin 
depth  in  the  soil  (a  «  6).  The  results  are  presented  for  an  incident 
exponential  pulse  that  propagates  as  a  plane  wave  that  is  uniform 
(constant  magnitude)  over  several  soil  skin  depths  along  the  cable.  The 
skin  depth  6  in  the  soil  is  plotted  in  Figure  2  as  a  function  of  soil 
conductivity  and  frequency.  The  decay  time  constant  t  of  the  exponential 
pulse  is  large  compared  to  the  soil  time  constant  =  eQ/cr.  The  small 
resistance  of  the  cable  is  neglected,  and  it  is  assumed  that  the  depth 
of  burial  is  small  compared  to  the  soil  skin  depth  at  the  highest 
frequencies  of  interest. 

2 .  Response  Characteristics 

The  total  current  induced  far  from  the  ends  of  a  long  cable  by 
-t/T 

an  incident  field  E  e  is  given  by 

o 


I  (OJ)  ~ 


I 


_ 1 _ 

O  /  jlUT  (jU)  +  1/t  ) 


(14) 
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in  the  frequency  domain,  or 


i  (t) 


I 


-t/T 

e 


o 


_2_ 


2 

u 


e 


du 


(15) 


it 

ible 


t 

rhe 

tial 

11 

h 


by 


(14) 


in  the  time  domain,  where 

g 

I  ^10  /t  t  E  D(t|r,cp)  (amperes) 

o  e  o 

e  -12 

o  8.85  X  10  , 

t  =  —  =  -  =  Time  constant  of  soil  (seconds) 

e  a  a 

t  =  Decay  time  constant  of  incident  pulse  (seconds) 

E  =  Peak  electric-field  strength  of  incident  pulse  (v/m) 

o 

D(i|r,cp)  =  cos  cp  for  vertically  polarized  wave 

=  sin  \jj  sin  cp  for  horizontally  polarized  wave. 

The  waveform  for  i(t)  and  the  incident  exponential  pulse  are  shown  in 
Figure  3  as  a  function  of  time  in  incident-pulse  decay  time  constant. 

The  peak  current  induced  in  the  cable  is 


pk 


=  I 


-t/T 


xj; 


V  t/T  2 
u 

e  du 


=  0.61  I 

max  o 


(16) 


and  the  peak  current  occurs  at 

t  =  0.85  T  .  (17) 

pk 

The  variation  of  the  peak  current  with  azimuth  angle  of  incidence  cp  and 
elevation  angle  of  incidence  i|i  is  shown  in  Figure  4.  The  magnitude  of 
the  peak  current  for  maximum  coupling  [D(\|r,cp)  =  l]  is  plotted  in  Figure 
5  as  a  function  of  soil  conductivity  for  various  incident-field  time 
constants.  The  plots  for  constant  t  are  truncated  in  the  lower  left-hand 
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TIME  — t  It 

FIGURE  3  WAVEFORMS  OF  INCIDENT  EXPONENTIAL  PULSE  AND  CABLE  CURRENTS 


SA-2192-4 


part  of  the  graph  where  t  =  t  .  The  responses  given  in  this  section 

are  valid  only  if  r  >  t  .  Where  the  approximations  used  here  are  valid, 

e 

the  peak  current  is  inversely  proportional  to  the  square  root  of  the  soil 

conductivity  (i. e. ,  sfo  I  =  constant). 

o 

Example:  Find  the  current  induced  in  a  buried  cable  by  a 

horizontally  polarized  exponential  pulse  of  10  kV/m  with  a 
decay  time  constant  of  lp,s  incident  from  an  elevation  angle  of 

10 


RELATIVE  MAGNITUDE  OF  PEAK  CURRENT  D(l p,(p) 
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FIGURE  4  VARIATION  OF  PEAK  CABLE  CURRENT  AS  AZIMUTH  (ip)  AND  ELEVATION  W 


PEAK  CURRENT,  lpEAK  —  kA/kV/m 


FIGURE  5  PEAK  CABLE  CURRENT  AS  A  FUNCTION  OF  SOIL  CONDUCTIVITY  AND 
INCIDENT  EXPONENTIAL  PULSE-DECAY  TIME  CONSTANT  r  =  U 
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30  degrees  and  an  azimuth  angle  of  70  degrees.  The  soil 
conductivity  is  10  mho/m. 

-3 

From  Figure  5  we  find  that  for  a  =  10  and  T  =  1  [is , 

_2 

the  peak  current  is  5.7  X  10  kA/kV/m  when  D(\|/,cp)  =  1.  There¬ 
fore,  for  10  kV/m,  the  peak  current  will  be  0.57  kA  when 
D(i|r,cp)  =  1.  From  Figure  4,  we  find  that  for  horizontal 
polarization  with  \|j  =  30°  and  cp  =  20°,  D(\|j,ep)  =  0.47. 

Therefore  the  peak  current  for  our  angle  of  incidence  will  be 
0.47  X  0.57  =  0.27  kA,  and  the  peak  will  occur  0.85  T  =  0.85  ^s 
after  the  beginning  of  the  current  pulse. 


3.  Theoretical  Basis 


-t/T 

The  incident  electric  field  E  e 

o 

Laplace)  transform 


has  the  Fourier  (or 


E(cju) 


E 

o 

jtu  +  1/T 


(18) 


The  resultant  z-component  of  this  field  at  the  surface  of  the  ground  for 
horizontal  polarization  is 


E 

z 


(1  +  R  )  E (cu)  sin 
h 


9 


E(o))  sin  i|r  sin  cp 


(a  »  cue  )  (19) 

o 


and  for  vertical  polarization  it  is 


E  =  (1  -  R  )  E(a>)  cos  cp  (a  »  cue  )  (20) 

z  v  o 


»  2  /—  / jOJ  E(cu)  cos  cp  (a  »  (JueQ) 


where  R  and  R  are  the  reflection  coefficients  for  oblique  incidence 
h  v 

on  the  air/ soil  interface. 
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\J2  6 


where  a  is  the  radius  of  the  cable  sheath,  y  =  1.781  ...  ,  and 

o 


1  +  J 


a  »  (Jue 


(24) 


To  obtain  the  final  expression  above,  use  was  made  of  the  fact  that 

log  (//26/y  a)  10  for  a  wide  range  of  soil  conductivities,  conductor 
o 

radii,  and  frequencies  (see  Figure  6). 


FIGURE  6  VARIATION  OF  LOG  (y/28/y0a)  WITH  FREQUENCY  AND  CONDUCTOR  RADIUS. 

(Note  that  a  factor  of  10  in  conductor  radius  produces  the  same  effect  as  a  factor 
of  100  in  soil  conductivity.) 


Using  this  expression  for  Z  in  Eq.  (22)  and  dropping  the 
exponential  term,  the  current  becomes 


rrE 


I(z,cju)  = 


jcup.  5 
o 


2tt/T~ 

_ e_ 

5jjl  /  ju) 
o 


E  D(i|t,cp) 
o 

jo)  +  1/t 


(25) 
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r 


J4) 


192-7 

S. 


25) 


Visualizing  the  buried  cable  as  a  long  (or  infinite)  lossy 
transmission  line,  the  current  far  from  the  ends  of  the  line  is 

I(z,u;)  =  P(z)e  YZ  +  Q(z)eYZ  (21) 


where 


E  z  z 

z  r  yz  -jk 


P(z)  =  J  e '  e 


2Z 

o 


dz 


Q(z)  = 


E  i  , 

Z  p  ^-yz  -jk 

97.  j 


e  e  z  dz 


o  z 


For  |  yz  |  »  1  and  yC&  -  z)  »  1  (i.e.,  for  points  several  soil  skin 
depths  from  either  end  of  the  cable),  the  current  is 


I  (z,cu) 


E 

z 

Z 


e 


-Jk 


/ 

z 


(22) 


The  exponential  term  gives  the  phase  of  the  incident  field,  which  is 
not  important  in  this  discussion.  The  impedance  Z  is  the  impedance  per 
unit  length  of  the  cable  and  soil  as  a  transmission  line.  This  impedance 
is  dominated  by  the  impedance  of  the  soil,  which  is  given  by 


Z  <=* 


2rraa 


(jya) 

(jya) 


Her  6 


+  j  log 


/26 


ya|  «  1 


5 

o 

TT 


(23) 
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where  T  =  e  / a ,  and  D(rt,cp)  is  the  directivity  function  given  in  Section 
e  o 

II-B-2  above.  The  inverse  transform  of  the  frequency-dependent  function 
in  Eq.  (25)  is 


-1 


/juj  (jou  +  1/T) 


=  /t 


-t/T 


_2_ 

0 


du 


(26) 


so  that 


i(z,t) 


/ T  T 
e 


E  D(i|i,cp) 
o 


-t/T 

e 


2 

u 

e 


du 


(27) 
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when  |jl  =  4tt  X  10 
o 

Useful  insight  into  the  physics  of  the  coupling  process  can 

be  obtained  by  using  the  second  expression  in  Eq.  (23)  for  the  impedance 

of  the  soil  and  noting  that  log  (/2  6/y  a)  »  tt/4  .  Putting  this  value 

o 

in  Eq.  (22)  and,  again,  dropping  the  exponential  phase  term,  the  current 
in  the  cable  is  seen  to  be 


2 

<jE  tt6 

Hz,w)  ~  Z  yf l  ■  <28) 

J  l0g  — 
o 

Note  that  aE  is  the  current  density  in  the  soil  (if  the  cable  were 
Z  2 

absent)  and  tt6  is  the  area  of  a  circle  one  skin-depth  in  radius  in  the 
soil.  Thus  the  total  current  induced  in  the  cable  is  approximately 
proportional  to  the  current  that  would  flow  in  a  circular  cylinder  of 
soil  one  skin  depth  in  radius  if  the  cable  were  not  present. 

For  cables  with  thin  insulation  over  the  shield,  the  impedance 
Z  per  unit  length  is 
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I 


z 


1 


5) 


O 


e 


t 


8) 


:e 


=  Z  +  jU)L 
is 

=  R  +  jUQ(L  +  L) 
s  s 


where,  for  6/a  »  1, 


R 

s 


L  +  L 

s 


JL 

2rr 


log 


y/26 

y  a 
o 


(29) 


(30) 


(31) 


A  plot  of  R/(d (L  +  L)  and  L  /L  in  Figure  7  illustrated  the 
s  s 

relative  magnitudes  of  the  resistance  and  reactance  of  the  impedance 

per  unit  length.  The  ratio  of  the  inductance  due  to  the  insulation  to 

that  due  to  the  soil  outside  the  insulation  in  Figure  4  shows  that 

most  of  inductance  is  attributable  to  the  soil  even  at  high  frequencies 

(6/a  small),  and  for  a  typical  insulating  jacket  thickness  of  the  order 

of  10  percent  of  the  cable  radius,  over  90%  of  the  inductance  is 

attributable  to  the  soil.  For  cables  a  few  centimeters  in  diameter  and 

soil  skin  depths  of  10  meters  or  more,  the  resistance  is  less  than 

one-fifth  the  inductive  reactance.  Therefore,  for  many  practical  cases 

we  may  use  the  approximation  R  «  <ju(L  +  L  )  to  simplify  the  solutions 

s 

to  the  buried-cable  analysis.  The  propagation  factor  for  the  cable/earth 
transmission  line  is 


Y 


+  jco(L 

s 


(32) 


Because  C(L  +  L)  »  CL,  the  propagation  factor  is  large  compared  to 
s 

k  cos  cp  cos  ^ ,  and  the  results  obtained  for  the  bare  cable,  assuming 
I Y |  »  k,  are,  to  a  good  approximation,  also  applicable  to  a  cable  with 
a  thin  insulating  jacket. 
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SKIN  DEPTH,  5/a  —  conductor  radii 
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FIGURE  7  VARIATION  OF  L/Ls  AND  R/co(L  +  Ls)  WITH  5/a  FOR  AN  INSULATED 
BURIED  CABLE.  (L  is  the  inductance  of  the  insulation  gap,  Ls  is  the 
inductance  of  the  soil,  and  R  is  the  soil  resistance,  all  on  a  per-unit-length 
basis.) 


C.  Current  Near  the  End  of  a  Long  Bare  Cable 
1.  Class  of  Problems 

The  response  characteristic  presented  here  are  for  points 
near  the  end  of  a  long  (semi-infinite)  bare  cable.  The  bulk  current  in 
the  cable  is  based  on  the  assumption  that  the  cable  sheath  is  in  contact 
with  the  soil  and  that  the  end  of  the  cable  is  terminated  in  a  very  low 
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0" 

-8 


impedance  (short  circuit)  or  a  very  high  impedance  (open  circuit) 
compared  to  the  characteristic  impedance.  The  short-circuit  case 
might  be  representative  of  a  cable  that  is  terminated  in  a  large  counter¬ 
poise  or  similar  low-impedance  structure.  The  open-circuit  case  is 
representative  of  a  cable  shield  that  is  dead-ended  or  insulated  at  the 
end  of  the  cable.  The  assumptions  that  a  >  cue  and  a  <  6  used  for  the 
long  bare  cable  above  apply  to  these  cases  as  well.  The  results  are 
presented  for  an  exponential  pulse  with  decay  time  constant  T  that  is 
large  compared  to  the  soil  time  constant  Tg  =  eQ/o.  The  electromagnetic 
pulse  is  incident  from  a  direction  defined  by  an  elevation  angle  i|f 
measured  from  horizontal  and  an  azimuth  angle  cp  measured  from  the  axis 
of  the  cable  (see  Figure  1). 

2.  Response  Characteristics 

The  total  current  induced  near  the  end  of  a  long  cable  that  is 
short-circuited  is  identical  to  the  current  far  from  the  ends  and  is 
given  by  Eqs.  (14)  and  (15),  and  has  the  waveform  shown  in  Figure  3. 

The  total  current  induced  near  the  end  of  a  long  cable  that 

-t/T 

is  open-circuited  at  the  end  by  an  incident  exponential  pulse  EQe 


is  given  by 


-  /  ju>/T  z/c 
e 

1  -  e 


l(z,0))  «  Io  /  jurr  (jcu  +  j/t) 


(33) 


in  the  frequency  domain,  and  by 


i  (z,  t)  <=» 


-t/T  i/W,  - 

°e  s*j  0  ' 


/  21  2 
»/u  I  e"  du 


(34) 


in  the  time  domain,  where 
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I  =10  /t  t  E  D(t|i.cp)  (amperes) 
o  e  o 


T  =  —  =  Time  constant  of  soil  (seconds) 
e  a 


t  =  Decay  time  constant  of  incident  pulse  (seconds) 

E  =  Peak  electric-field  strength  of  incident  pulse  (V/m) 
o 

D(i|r,cp)  =  cos  cp  for  horizontal  polarization 

=  sin  sin  cp  for  vertical  polarization 
z  =  Distance  from  end  of  cable  (meters) 


c  = 


1/v/TT  e  =  Speed  of  light  in  free  space 
o  o 


/  \  2 

\c/  4tt 


This  current  is  zero  (very  small)  at  the  end  (z  —  0)  and  approaches  the 

value  given  by  Ea.  (15)  at  large  distances  from  the  end.  Plots  of  the 

waveform  at  distances  of  3.16,  10,  and  31.6  meters  from  the  end  of  the 

-2 

cable  are  shown  in  Figure  8  for  a  soil  conductivity  of  10  mho/m. 

1  2  2  2 

Since  the  waveform  is  determined  by  (z/c)  =  az  / e  c  ,  the  waveforms 

Te  O 

shown  in  Figure  8  will  represent  different  distances  from  the  end  of 

-2 

the  cable  if  the  soil  conductivity  differs  from  10  mho/m.  The 

magnitude  of  the  current  remains  inversely  proportional  to  the  square 

root  of  the  soil  conductivity  (/ ~a  I  ^  constant). 

o 


3. 


Theoretical  Basis 


The  component  of  the  electric  field  along  the  cable  is  given 
by  Eqs .  (19)  and  (20).  This  field  induces  a  current  in  the  semi-infinite 

cable  extending  from  z=0toz=4-*°°,  which  is  given  by  Eq.  (8), 
where 
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1  when  shield  is  open-circuited 


=  -1  when  shield  is  short-circuited 
k7  =  k  cos  i|r  cos  cp. 

The  current  near  the  end  at  z  =  0  then  becomes 


I  (z,co) 


E  r  t,' 

z  -jk  z 

—  e 

Z  _ 


for  the  open-circuited  cable,  and 


(35) 


I(z,u>)  =  — 


E 

z  -Jk  z 
e 
Z 


(36) 


for  the  short-circuited  cable.  The  term  exp  (-k7z)  =  exp  (-jkz  cos  cp  cos 
gives  the  phase  of  the  incident  wave  at  the  point  z  relative  to  the  phase 
at  the  end  z  =  0.  Since  k7  «  |y|,  the  currents  become 


I  (z,w) 


E 

z 

Z 


(x  -  e--) 


(open  circuit)  (37) 


I(z,cd)  » 


(short  circuit)  (38) 


when  the  phase  is  referred  to  the  phase  of  the  locally  incident  field. 

The  short-circuit  current  is  thus  the  same  as  the  current  a  great 

distance  from  the  end  [see  Eq.  (25)].  The  current  near  the  open- 

circuited  end  of  the  cable  is  this  current  reduced  by  (1  -  e  ^Z),  where 

y  ^  y/jiup-  a. 
o 

For  the  exponential  incident  pulse  given  by  Eq.  (18)  and  the 
impedance  per  unit  length  given  by  Eq.  (23),  the  current  near  an  open 
circuit  is  given  by 
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I  (z,io) 


2 it/  t  t  E  D(f,cp) 

e  o  1 


jlJU/T  z/c 
e 


/  JU)T  <  jCJL>  +  1/t) 


where  /u,  a  is  written  as  1/c  /t  .  The  first  term  in  the  frequency- 
o  e 

dependent  part  of  Eq.  (39)  is  similar  to  Eq.  (14)  and  has  an  inverse 
transform  similar  to  Eq.  (15).  The  second  term  can  be  written  as  a 
convolution  of  the  exponential  pulse  with  the  impulse  response  of  the 
function  containing  exp  /  ju>  /  /"jiiu  obtained  from 


„  -/  jU)/T  z/c 

(J0  ■“  1  e  g 

/  jO)T 


exp  -  (z/c)' 


1 

T  t  / 

e  J  /  1 


v% « 


SO  .  (40) 


The  current  waveform  is  then  given  by 


Ife.t)  =  Ioe-t/J  .  e-P/“2)/  du  (41) 


where 


I  =10  /n  E  D(ilt ,cp) 
o  e  o 


(£)2  — 

\  c  /  4t 
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Ill  CURRENTS  AND  VOLTAGES  INDUCED  ON  INTERNAL  CONDUCTORS 


A.  General 

The  internal  currents  and  voltage  in  this  section  are  presented  for 

exponential  currents  in  the  shield  of  a  tubular  shielded  cable.  In 

Section  III-B,  the  internal  currents  and  voltages  are  given  for  cables 

that  are  short  compared  to  the  shortest  wavelength  penetrating  the  shield. 

In  Section  III-C,  these  results  are  extended  to  cables  that  are  long 

compared  to  the  shortest  penetrating  wavelength,  neglecting  attenuation 

of  signals  propagating  on  the  core.  In  both  sections,  the  results  are 

presented  for  shield-current-decay  time  constants  t  that  are  short, 

intermediate,  and  long  compared  to  the  shield  diffusion  time  constant 

t  .  In  Section  III-D,  some  results  are  presented  for  the  case  in  which 
s 

shield  currents  induced  in  a  buried  cable,  determined  from  a  wave 
incident  on  the  surface  of  the  ground,  are  used  to  compute  the  internal 
currents. 

B.  Electrically  Short  Cable  with  Tubular  Shield  (Cable  Current  Given) 

1 .  Class  of  Problems 

The  type  of  shield  considered  here  is  a  tubular  shield  with  no 
holes  and  no  seams  that  are  not  continuously  welded,  brazed,  or  soldered. 
The  results  are  not  applicable  to  braided-wire  or  tape-wound  shields. 

The  cables  are  electrically  short  at  the  shortest  wavelength  penetrating 
the  shield  (see  Figure  9).  The  results  are  given  for  zero-rise-time 
exponential  pulses;  however,  these  results  are  applicable  to  finite-rise¬ 
time  pulses  if  the  pulse  rise  time  is  short  compared  to  the  diffusion 
time  (see  Figure  10).  The  open-circuit  voltage  between  the  core  and  the 
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FIGURE  10  THE  DIFFUSION  CONSTANT  ts  AS  A  FUNCTION  OF  WALL  THICKNESS  FOR 
COMMON  SHIELDING  MATERIALS 
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shield  at  the  ends  of  the  cable  is  presented  for  a  given  exponential 
pulse  of  current  flowing  in  the  cable  shield.  Also  presented  is  the 
current  in  the  core  when  the  terminating  impedance  between  the  core  and 
shield  at  each  end  is  the  common-mode  characteristic  impedance  of  the  core- 
to-shield  transmission  line. 


2.  Response  Characteristics 


For  a  cable  of  total  length  &  whose  shield  has  a  dc  resistance 

—  t/T 

per  unit  length  R  and  a  current  I  e  flowing  in  it,  the  open-circuit 

o  so 

voltage  developed  between  the  core  and  the  shield  at  the  end  of  the  cable 
is 


V  (uo) 


I  R  i 
so  o 


2 


/jurT 

s 


sinh  /jurr 

s 


(42) 


where 

R  =  (2iraaT)  ^  =  dc  resistance  per  meter  of  the  shield 
o 

2 

T  =  {.lctT  =  Diffusion  time  constant  for  the  shield  (see  Figure  10)) 
s 

r  =  Decay  time  constant  of  the  exponential  pulse  of  shield  current. 

The  quantities  a,  cr,  \±,  and  T  are  the  radius,  conductivity,  permeability, 
and  wall  thickness  of  the  shield,  respectively.  The  voltage  waveform  is 


V(t)  =  V  f(t) 
o 


where 


f(t) 


/jurT 


^ju>  +  sinh  /  juST 


and  V 


I  R  l 
o  o 

2 
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The  function  F(t)  is  plotted  in  Figure  11  for  t  =  Tg  and  for  T  »  Tg. 

For  t  «  t  the  function  f^ (t)  (t^/t)  is  plotted.  The  coefficient 
(per  kilometer)  is  plotted  in  Figure  12  as  a  function  of  shield  thickness 
and  diameter  for  copper,  aluminum,  steel,  and  lead  shield.  The  peak 
open-circuit  voltages  are 


v(t) 


peak 


=  0.77  V 

o 


(T  »  T  ) 
S 

(T  «  T  ) 
S 

(T  =  T  ) 

S 


(43) 


if  the  core  is  open-circuited  (relative  to  the  shield)  at  both  ends 
The  10-to-90-percent  rise  time  for  the  voltage  is 


t  =  0.236  T  (T  >;>  Ts>  <44> 

10-90  s  fa 

=  0.038  T  (T  «  T  ) 

s  s 


=  0.15  T 

s 


(T  =  T  ) 

s 


Plots  of  the  asymptotic  values  of  the  peak  voltage  and  rise  time  against 

t/t  are  shown  in  Figure  13  as  solid  lines,  and  an  estimate  of  their 
s 

behavior  between  asymptotes  is  shown  as  a  dashed  curve. 

The  voltage  has  the  same  magnitude  and  shape  at  both  ends  of 
the  cable;  however,  the  polarity  at  one  end  is  opposite  to  the  polarity 
at  the  other  end.  If  the  core  is  shorted  to  the  shield  at  one  end,  the 
open-circuit  voltage  at  the  opposite  end  will  have  the  same  waveform  but 
it  will  double  in  magnitude. 

The  current  through  matched  terminations  at  the  ends  of  the 
cable  will  be 


Kuo)  =  V(uo)/Z 

o 


(45) 
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o  0.1  0.2  0.3  0.4  0.5 

TIME  — t/Ts 

(a)  IMPULSE  RESPONSE  [Tl$05(t)] 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

TIME  —  t/T$ 

(b)  STEP  FUNCTION  RESPONSE  [ISQu(t)) 


FIGURE  11  CORE  CURRENT  WAVEFORMS  PRODUCED  BY  AN  EXPONENTIAL  PULSE 
L„e_t/T  OF  CURRENT  IN  THE  SHIELD  (normalized  to  V0/Z0  =  IsoR0^2Zo' 
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FIGURE  12  OPEN-CIRCUIT  VOLTAGE  V0  =  IS0R02/2  BETWEEN  CORE  AND  SHIELD  FOR  COMMON  SHIELDING 
MATERIALS,  THICKNESSES,  AND  DIAMETERS  (for  Iso  =  1  ampere  and  8  =  1  kilometer) 
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or 


i(t)  =  v(t)/Z 

o 

where  Z  is  the  common-mode  characteristic  impedance  of  the  core-to-shield 
transmission  line.  This  current  will  be  the  same  magnitude,  waveshape, 
and  polarity  throughout  the  length  of  the  cable  for  the  electrically 
short  cables  considered  here.  Thus,  all  of  the  data  for  the  open-circuit 
voltage  can  be  used  to  obtain  core  current  if  the  characteristic  impedance 

Z  can  be  obtained, 
o 

Example:  A  lead-sheathed  cable  2  inches  in  diameter  with 
a  100-mil-thick  sheath  is  two  miles  (3.2  km)  long  and  is 
subjected  to  an  exponential  sheath  current  of  1000  amperes 
whose  decay  time  constant  is  1  [is .  Find  the  peak  open 
circuit  voltage  developed  at  the  ends  of  the  cable. 

From  Figure  9,  T  =  6  for  the  lead  sheath  when  f  =  9  kHz 

or  X  =  33  km,  so  the  3.2-km  segment  is  electrically  short 

for  the  penetrating  frequencies.  From  Figure  10,  the 

diffusion  time  constant  t  =  3,5  X  10  seconds,  which  is 

s  -2 

much  larger  than  t,  and  t/t  =  2.86  X  10  .  Therefore,  the 

s 

peak  voltage  is  approximately  (T/Tg)  Vq.  From  Figure  12, 

is  0.55  volts/km  for  1  ampere,  or  V  =  0,55  X  1000  X  3.2 

o 

^  1760  volts,  and  from  Figure  13  the  peak  voltage  is  about 
0.16  V  p-  280  volts.  The  10-to-90-percent  rise  time  from 
Figure  13  is  0.038  t  =  1.33  X  10  seconds.  The  waveform 
would  be  that  shown  in  Figure  11(a). 

3.  Theoretical  Basis 

The  thin-walled  tubular  shield  of  radius  a,  wall  thickness  T, 

conductivity  a,  and  permeability  u  =  u,  ll  has  a  transfer  impedance  given 

r  o 
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by 


1 


/jujicr  T 
sinh  y  jajjjL.O'  T 


Z  = 

T  2itaaT 


/  Jourr 


=  R 


o  sinh  /  j(JDT 


(46) 


where  R  =  (2itaaT)  1  is  the  dc  resistance  per  unit  length  of  the  shield 
°  2 

and  Tg  =  haT  is  a  time  constant  for  diffusion  through  the  shield.  For 

cables  that  are  short  compared  to  the  shortest  wavelength  that  penetrates 

the  shield,  the  current  induced  in  internal  conductors  that  are  terminated 

in  their  characteristic  impedance  Z  at  both  ends  of  the  cable  by  a 

o 

current  I  flowing  in  the  shield  is 
s 


1(0)) 


IRA  /jorr 
s  o  _ s 

2Z  sinh  /  jorr 
o  s 


(47) 


where  A  is  the  length  of  the  cable. 

The  length  A  for  which  this  estimate  is  valid  can  be  ascertained 

from  Figure  10,  where  the  diffusion  time  constant  t  is  plotted  as  a 

s 

function  of  wall  thickness  for  various  shielding  metals.  When  t  >  A/c, 

s 

the  diffusion  time  is  greater  than  the  transit  time  along  the  length  of 
the  cable,  and  the  entire  cable  adjusts  to  the  penetrating  signal  as  it 
is  developed.  When  the  diffusion  time  is  less  than  the  transit  time, 
the  propagation  times  to  various  parts  of  the  cable  become  significant 
and  transmission-line  theory  must  be  used.  An  alternative  method  of 
determining  the  validity  of  the  short-line  approximation,  which  is 
applicable  to  the  frequency  domain,  makes  use  of  Figure  9,  which  is  a 
plot  of  the  frequency  at  which  the  skin  depth  in  the  shield  is  equal  to 
the  wall  thickness.  When  T/6  >  lf  the  attenuation  by  the  shield 
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3S 

ted 


) 


.ned 


f 

t 


is  large,  whereas,  when  T/6  <  1,  the  attenuation  is  small.  Therefore, 
the  wavelength  at  which  T/6  =  1  is  the  shortest  wavelength  that  freely 
diffuses  through  the  shield.  The  free-space  wavelength  associated 
with  the  frequency  at  which  T/6  =  1  is  indicated  on  a  scale  along  the 
top  of  Figure  9. 

When  the  shield  current  is  an  exponential  pulse  of  the  form 

_+-/t 

I  e  ,  whose  transform  is 
so 


I  (cu) 
so 


I 

so 

jU)  +  1/t 


the  current  in  the  internal  conductors  is 


I(uo)  = 


I  R  i 

SO  o 

2Z 


/jUJT 


^juu  +  ~  j  sinh  /durrt 


(48) 


Note  that  if  t  »  T  ,  the  significant  spectrum  of  I(oj)  is  determined  by 

s 

jivr  ,  and  ju)  +  1/t  «  jo)  in  the  denominator.  The  response  of  the  internal 
s 

conductors  is  then  the  step- function  response  of  the  shield.  If,  on 

the  other  hand,  T  «  t  ,  the  significant  spectrum  of  1(a))  is  again 

1 

determined  by  jurr  ,  but  now  juo  +  1/t  ^  —  in  the  denominator.  The  response 
s  T 

of  the  internal  conductors  is  then  the  impulse  response  of  the  shield 

if  the  impulse  is  represented  by  Tl  6(t),  where  6(t)  is  the  unit 

so 

impulse. 

The  response  waveform  of  the  internal  current  for  the  short 

pulse  (t  «  t  )  is  thus 
s 


i(t)  = 


n=l 


(2n-l)  t 


2t 


-  1 


exp 


-(2n~l) 


T 

2  _ s 

4t 


(T  «  T  ) 

s 


(49) 
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where  I  „ 
i 


I  Rl 
so  o 

2Z 


o 


f.(t) 

0 


and  for  the  long  pulse  (t  »  t  )  it  is 

s 


i(t) 


1  u  „  A  V  •  r  t  n 

e-p[w^] 

o  '  '  n=l  L  J 


(T  »  T  )  (50) 

s 


=  I,  f  (t) 
a  u 


The  functions  f  (t)  and  f  (t)  are  plotted  in  Figure  11,  as  well  as  the 
6  u 

response  for  an  exponential  pulse  for  which  t  =  T  . 

s 

The  voltage  developed  between  the  internal  conductors  and  the 
shield  at  the  ends  of  the  cable  when  both  ends  are  open-circuited  is  Z 

o 

i(t) ,  or 


(t  «  t  )  (51) 

s 

(T  »  T  )  (52) 

s 

C.  Correction  Factor  for  Long  Cables 
1.  Class  of  Problems 

The  results  given  below  permit  the  peak  core  current  and  time- 
to-peak  to  be  determined  for  cables  of  any  length.  The  results  are 
applicable  to  tubular  shields  in  which  an  exponential  current  pulse 
propagates  in  one  direction  only  (no  reflections).  The  core-to-shield 
transmission  line  is  assumed  to  be  lossless  (no  attenuation)  and  terminated 
in  its  characteristic  impedance  at  its  ends. 


v(t)  = 


I  R  i  /  v 
SO  O  /  T  \ 

2  Iv 


f.(t) 

0 


I  R  i 
so  o 


f  (t) 
u 
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2. 


Response  Characteristics 


0) 


ie 

7 

a 

o 


51) 


52) 


me- 


The  peak  current  through  matched  terminations^ 
shielded  cable  of  any  length  l%  when  a  current  1^  e 
shield,  is,  for  t  «  t  , 


at  the  end  of  a 
1  flows  in  the 


i 


peak 


I  R  i 

50  °-  JL  F 

2Z  T  H 

o  s 


(T  «  T  )  (53) 

s 


where 

F  —  5,9 

l 


T 

S 

~  t 

O 

t  =  —  C/e"  ±  k//k) 
o  c  r 


e  =  Dielectric  constant  of  core-to-shield  insulation 


t  =  Shield  diffusion  time  constant  (see  Figure  10) 


k'/k  =  1,  if  shield  current  is  propagating  with  velocity  c 


=  cos  ep  cos  if  shield  current  is  induced  in  a  buried 
cable  by  a  wave  incident  from  (ty,cp)  (see  Figure  !)• 


The  positive  sign  is  used  when  the  direction  of  propagation  of  the  shield 
current  is  away  from  the  terminals,  and  the  negative  sign  is  used  when 
the  direction  of  propagation  is  toward  the  terminals  at  which 
to  be  determined,  A  plot  of  F  as  a  function  of  t  /t  is  shown  in  Figure 


.d 

ainated 
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LENGTH  CORRECTION  FACTOR, 


FIGURE  14  CORRECTION  FACTOR  FOR  THE  PEAK  VOLTAGE  AND  CURRENT  INDUCED 
IN  LONG  SHIELDED  CABLES  (k'  is  the  propagation  factor  for  the  shield  current, 
k  =  co\ZfX0e0  ,  and  er  is  the  dielectric  constant  of  the  core-to-shield  insulation) 


The  time  to  reach  the  peak  currents  is 


t  ,  «  0,09  t 
peak  s 


(55) 


^  1*05  t 

o 


0.5  T 

s 
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rent, 


55) 


The  waveform  of  the  current  is  approximately  that  of  Figure  11(a)  for 

t  /T  <o.l,  and  like  the  leading  edge  of  Figure  11(b)  for  >0.5 

o  s 

Typical  waveforms  are  shown  in  Figure  15(a). 

The  peak  voltage  across  the  matched  terminations  is  1peak 
z  if  both  ends  of  the  cable  are  terminated  in  the  characteristic 
impedance  Z  .  The  peak  open-circuit  voltage  at  one  end  when  the  other 
end  is  terminated  in  its  characteristic  impedance  is  2  ipeak  Zq. 

When  t  »  t  ,  the  peak  current  in  the  core  is 


I  R  H 
so  o 

1peak  ~  2Z 

o 


(T  »  t  )  (56) 

o 


and  the  time  to  reach  the  peak  (shoulder  of  the  current  step)  is 


t  fa  0.5  T 
peak  s 


S3  t 

O 


(t  «  T  ) 

O  S 

(57) 

(t  <5<  T  ) 
o  s 


Waveforms  for  the  core  current  are  shown  in  Figure  15  for  T  <<:  Tg> 

t  »  t  ,  and  t  =  t  ,  in  terms  oft  =  l/v.  The  transit  time  A/v  is 
o  S  O 

expressed  in  terms  of  the  10-to-90-percent  rise  time  of  the  step- 

function  response  of  a  short  cable.  The  value  of  the  10-to-90-percent 

rise  time  is  t  =  0.236  T  .  A  plot  of  t  for  polyethylene  insulated 
r  s  o 

cable  (e  =2.3)  and  various  values  of  kVk  =  cos  cp  cos  is  shown  in 
r 

Figure  16.  Note  that  kVk  =  1  gives  the  condition  for  the  shield  current 
propagating  with  velocity  c,  regardless  of  whether  the  current  is  induced 
by  an  incident  wave  or  by  some  other  means. 


3.  Theoretical  Basis 

A  cable  with  a  tubular  shield  is  assumed  to  extend  from 
z  =  0  to  z  =  A,  and  its  core  is  terminated  in  its  characteristic  impedance 
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CABLE  LENGTH,  £  —  km 

SA-2192-17 


FIGURE  16  DIFFERENTIAL  TRANSIT  TIME  tQ  AS  A  FUNCTION  OF  CABLE  LENGTH 
AND  PROPAGATION  FACTOR  k'  FOR  THE  SHIELD  CURRENT  FOR 
POLYETHYLENE-INSULATED  CABLE  CORE 
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at  both  ends.  The  shield  current  is 


T  ±jk'z 

I  e  , 

so  _  ±jk  z 

~  T  I  0 

jo)  +  1/t  so 


(T  «  t  )  (58) 

s 


and  the  core  current  through  the  termination  at  z  =  0  is,  from  Eqs .  (8) 
through  (13), 


KO,  u))  =  Q(0) 


I  R 
so  o 

2Z 


o 


t  y  jojr 

_ s 

sinh  /jour 

s 


-jd)t 

O 


jUJt  /l 

O 


(59) 


where  tQ  =  —  ±  —  j  ,  is  the  dielectric  constant  of  the  insulation 

between  the  core  and  the  shield,  and  k  =  u)/c.  The  current  waveform  is 
thus 


I  R  JL 

i(0, t)  =  ~ °°  —  [f  (t)  -  f  (t-t  )  u(t-t  )] 

ZZi  t  u  u  o  o 

o  o 


I  Rl 
so  o 


2Z 

o 


—  F(t) 

T 


S 


(60) 


where 

T 

F(t)  =  [f  (t)  -  f  (t-t  )  u(t-t  )]  .  (61) 

t  u  u  o  o 

o 

The  function  f^(t)  is  plotted  in  Figure  11(b)*  From  Figure  11(b)  it  is 

apparent  that  for  0.05  ^  t/T  ^  0.5,  the  maximum  value  of  F(t)  will 

s 

occur  when  t/r^  ^  1.05  t^/x^,  and  the  value  of  this  maximum  will  be 

t^)/(t^/T^).  When  t^  «  t  .  the  value  of  F(t)  can  be  obtained 
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(60) 


from 


the  short-cable  form  of  Eq.  (59),  for  which 


>8) 


1  - 


juu  t  /A 
o 


(62) 


(8) 


and 


IRA 

v  SO  o  T 

i(0,t)  = - f 


2Z 


t  6  (t) 
s 


(63) 


59) 

ation 

is 


The  function  f  (t)  is  plotted  in  Figure  11(a)  and  has  a  maximum  value  of 
6 

5.9  at  t  =  0.09  t  .  Therefore,  the  value  of  F(t)  that  produces  the  peak 
s 

value  of  i(0,t)  is 


F,  «  5.9 
A 


<  0.1 


(64) 


f  (1 . 05 t  ) 
u  o 


—  >  0.05 
T 

S 


(61) 


ft 


—  >  0.5 

T 

S 


If  the  shield  current  is  propagating  with  the  speed  of  light, 


t  is 


k  =  k  and 


ned 


t  =  —  (flT  ±1) 

o  c  r 


(65) 
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A  more  interesting  case  for  buried  cables,  however,  is  the  case  in  whi 
the  shield  current  is  generated  by  the  local  incident  field.  In  this 
case,  k/  =  k  cos  if  cos  cd,  and 

l 

t  =  —  (/e~  ±  cos  if  cos  cp) 
o  c  r 

When  the  pulse  time  constant  t  is  large  compared  to  the  diffe 

ential  transit  time  t  ,  the  response  to  a  step  function,  which  is  the 

o 

integral  of  the  impulse  response,  is  obtained.  Thus 


1  I  R  l  t 

S  i(0,t)  =  2Z  ~  j  F(t)dt  <T  »  t  >  (6< 

1  o  0 

I 

> 

} 

}  and  the  step-function  of  shield  current  produces  a  step  function  of  in¬ 

ternal  current.  The  rise  time  of  the  internal  current  is  not  zero, 

however,  as  can  be  seen  in  Figure  15(b).  If  t  is  small,  the  approxims 

o 

of  Eq.  (62)  leads  to 


I  R  a 

i(0,t)  w  f  (t)  (t  «  T  )  (67 

2Z  u  os 

o 

while  if  t  is  not  small  compared  to  t  ,  the  integration  of  Eq.  (62)  le 
o  s 

to 


i (0, t) 


I  R  jfc 
so  o  t 

2Z  t 
o  o 


(t  ^  t  ) 
o 


I  Rl 
so  o 


(t  >  t  ) 
o 


(68 


Thus,  in  either  case  the  peak  current  is  I  R  1/2Z  ,  but  the  time  to 

so  o  o 

reach  the  peak  is  approximately  0.5  t  when  t  «  t  ,  and  is  approximate 

sos 
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t  when  t  is  not  small  compared  to  T  . 
o  °  S 


D>  internal  Voltage  from  Incident  Field 
1 <  Class  of  Problems 

The  results  presented  below  give  the  voltage  between  the  core 
and  the  shield  induced  in  a  buried  cable  by  an  exponential  plane-wave 
pulse  incident  on  the  surface  of  the  ground.  The  results  apply  to  cables 
that  are  short  compared  to  the  shortest  wavelength  penetrating  the  shield 
(see  Figure  9).  The  cables  are  assumed  to  have  at  least  one  tubular  shield 
(see  Section  III-E).  The  results  are  obtained  from  the  convolution  of 
the  cable  current  obtained  in  Section  III-B  with  the  impulse  response 
of  the  tubular  shield  obtained  in  Section  III-B. 


2.  Response  Characteristics 

The  open-circuit  voltage  induced  between  the  core  and  the  shield 

of  a  buried,  tubular-shielded  cable  by  an  incident  exponential  pulse 

E  e  is  given  by 
o 


_ 1 _ 

(ju)  +  1/t)  sinh  /  jurr 

s 


(69) 


where 


R  i 

V  =  -2-  [see  Eq.  (42)] 

o  z 

6 

I  =10  /t  t  E  D('|f,cp)  [see  Eq.  (15)] 
o  e  o 

2 

T  =  (j,<rT  is  the  diffusion  constant  for  the 
s 

tubular  shield. 
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The  open-circuit  voltage  waveform  is 


This  waveform,  normalized  to  is  plotted  in  Figure  17  for  several 

values  of  t  /t  . 


3.  Theoretical  Basis 

The  current  induced  in  a  buried  cable  by  an  incident  exponential 

field  E  e  is  given  by  Eq.  (14)  as 

o 


_ 1 _ 

I(u))  =  ^/Surf  (jorr  1/t) 


(71) 


and  the  open-circuit  voltage  developed  between  the  core  and  the  shield 
is  given  by  Eq.  (47)  as 


V((u)  =  Vo  sinh  /  jucrr 

s 


(72) 


for  an  impulse  of  shield  current.  The  open-circuit  voltage  induced  by  a 
shield  current  I(u))  is  therefore 


V(u>) 


_ _L _ 

(jtu+l/r  )  sinh  /  jurT 


(73) 


and  the  open-circuit  voltage  waveform  is 


i: 


1 

rite 


47 


v(t)  =  V  I 
o  o 


s  -t/r 
—  e 
nr 


(t  /r)y 
s 
e 


E 


(74) 

The  function  v(t)  normalized  to  V  I  /f rc  is  plotted  in  Figure  17  for 

o  o 

several  values  of  t  /t.  The  current  i(t)  through  matched  terminations  on 

5 

the  core  is  v(t)/Z  . 

o 


E.  Multiple  Shields 

1.  Class  of  Problem 

The  results  presented  below  apply  to  multiple-shield  cables 
having  at  least  one  tubular  shield.  The  remaining  shields  may  be  braided 
wire,  tape-wound,  or  any  other  leaky  shield  whose  transfer  impedance  is 
relatively  independent  of  frequency  in  the  range  of  frequencies  passed 
by  the  tubular  shield.  The  results  below  provide  rough  approximations 
for  equivalent  tubular  shields  having  the  same  shielding  properties  as 
the  multiple  shield.  These  equivalent  shield  characteristics  can  then  be 
used  with  the  results  of  Sections  III-B,  -C ,  or  -D  to  obtain  internal 
responses. 


2.  Equivalent  Shield  Characteristics 

A  multiple-shield  system  containing  one  tubular  shield  whose 

transfer  impedance  is  R  V  jojr  /sinh  /  jorr  and  other  shields  whose  combined 

01  s  s 

resistance  per  unit  length  is  R^  for  frequencies  between  0  and  l/2?t 
has  an  equivalent  transfer  impedance  given  by 


(R  «  R  )  (75) 

01  02 
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and 


l  - == 

02  sinh  /  jorr 


(R  »  R  )  .  (76) 

01  02 


Theoretical  Basis 


If  a  cable  shield  consists  of  several  layers  of  shielding 
material,  at  least  one  of  which  is  a  tubular  shield  with  no  holes  or 
seams,  the  transfer  impedance  of  the  composite  shield  can  be  written 

VT 

1  -  -  - 

_ 01  02  sinh  yT _  ^ 

T  ”  R  coth  vT  +  R  „  +  j(j)L 
01  02 


where  R  is  the  dc  resistance  per  unit  length  of  the  tubular  shield  of 
01 

thickness  T;  R  is  the  dc  resistance  per  unit  length  of  the  other 
02 

shields;  yT  =  /jurr  ;  and  L  is  the  inductance  per  unit  length  between 
s 

the  tubular  shield  and  the  other  shields.  It  is  assumed  that  the 

transfer  and  internal  impedances  are  relatively  independent  of  frequency 

when  yT  £  1,  and  it  is  assumed  that  o)L  is  relatively  insignificant 

compared  to  R  or  R  when  vT  £  1.  Then 
01  02 


R  R  yT 
01  02  ' 

jr  .  .  ' 

T  R  cosh  yt  +  r02  sinh  YT 


<“L  <<:  V  V 


Now,  when  R  „  »  R,. 
’  01  02 


yt 

2  R  1 

T  ~  02  cosh  yt 


and  since  cosh  yt  and  sinh  yT  behave  the  same  for  large  yt>  the  high- 
frequency  behavior  of  the  transfer  impedance  is  very  similar  to  that 
of  a  tubular  shield  whose  dc  resistance  is  R  and  whose  diffusion 

\JCi 
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characteristics  are  those  of  the  tubular  shield.  The  current  through 
matched  terminations  between  the  core  and  shield  is,  by  analogy  with 
Eq.  (48), 


IsoR02/j 

2Z  1  — 

°  (ju)  +  — )  cosh  /ju W 


(R  »  R  ) 
01  02 


where  Tg  is  determined  from  the  tubular  shield  parameters,  and  R  is 

02 

determined  from  the  other  shields.  The  current  waveforms  are  given  by 


i(t)  «  I 


T  1  T  7 

T  /it  T 
s  s  n=l 


(2n-l)  t 


-  -1  (-l)n  1eXp 


■<2“-1)2  if 


(T  «  T  ) 
S 


f  E 


exp  [-(2n-l)  t  /4t] 
s 


(T  »  T  )  (81) 

s 


I  R  i 
so  02 


w  ere  L ^  =  2Z  •  Comparing  these  expressions  with  Eqs.  (49)  and 

(50),  it  is  observed  that  the  series  above  alternates  because  of  the 
(  -»n_1 

(  1)  coefficient,  whereas  those  in  Eqs.  (49)  and  (50)  do  not 
alternate.  Furthermore,  it  has  been  found  that  the  first  term  of  the 
series  defines  the  waveform  quite  well  for  times  beyond  the  peak  of 
Figure  11(a)  and  beyond  the  shoulder  of  Figure  11(b).  Therefore  the 
alternating  terms  can  affect  only  the  late-time  behavior  of  the  wave¬ 
forms,  and  then  only  be  reducing  the  late-time  magnitude.  Therefore, 
the  waveforms  of  Figure  11  will  be  very  good  approximations  to  the 

waveforms  induced  through  the  composite  shield  with  R  »  R 

01  02* 
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Eq.  (77)  approaches 


R  vT 

Z  pa  - — -  (R  »R  ) 

T  sinh  yT  02  01 


(82) 


which  is  the  transfer  impedance  of  the  tubular  shield  alone.  In  this 
case,  the  other  shields  are  not  very  effective,  and  the  composite  shield 
behaves  as  though  it  were  a  single  tubular  shield.  The  results  of 
Section  III-B  then  apply  directly  to  this  area. 

The  quality  of  the  predictions  based  on  these  approximations 
depends  on  how  well  the  assumptions  are  met.  Generally,  the  assumptions 
are  of  such  a  nature  that  the  predicted  core  voltages  or  currents  will 
be  larger,  rather  than  smaller,  than  the  actual  voltages  and  currents 
induced  on  the  cable  core,  since  all  of  the  predictions  involve  reducing 
the  size  of  the  denominator  in  Eq.  (77),  which  increases  the  transfer 
impedance . 
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Appendix  A 


a 

c 

C 

D(\[f  ,cp) 

e 

E 

f 

G 

i 

I 

3 

k 

K 

i 

L 

P(z)  ) 
Q(z)  J 

r 

R 

R  ,  R 
v  1 

s 

t 

T 


v 

V 


LIST  OF  SYMBOLS 

Radius  of  cable  shield  (meters) 

Speed  of  light  (3  X  108  meters/second) 

Capacitance  (farads,  farads/meter) 

Directivity  factor  (dimensionless) 

2.718... 

Electric-field  strength  (volts/meter) 

Frequency  (hertz) 

Conductance  (mhos/meter) 

Instantaneous  current  (amperes) 

Monochromatic  current  (ampere-seconds) 

Imaginary  unit  J-l 

Free-space  phase  factor  (meter  ) 

Constant 

Cable  length  (meters) 

Inductance  (henries,  henries/meter) 

Source  functions  (amperes) 

Radial  distance  (meters) 

Resistance  (ohms,  ohms/meter) 

Reflection  factor  for  oblique  wave  incidence  (dimensionless) 
Complex  frequency  (second  ) 

Time  (seconds) 

Shield  thickness  (meters) 

Instantaneous  voltage  (volts) 

Monochromatic  voltage  (volt -seconds) 
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Y 
z 
Z 
Z 

o 

O' 

p 

Y 

Yo 

6 


e 

e 

o 

e 

r 

T1 

X 


u 

o 

Mr 

P 


a 

T 


9 

* 

a) 


Admittance  per  unit  length  (mhos/meter) 

Distance  along  axis  of  cable  (meters) 

Impedance  per  unit  length  (ohms/meter) 

Characteristic  impedance  Jz/Y  (ohms) 

Attenuation  constant  (nepers/meter) 

Phase  factor  k/e~  (meters  "S 

r  -1 
Complex  propagation  factor  (meters  ) 

1 , 781 . . ,  (Euler 1 s  constant) 

-  i 

Skin  depth  (irfj uo)  2  (meters) 

Permittivity  e  e  (farads/meter) 

r  °  -12 

Permittivity  of  free  space  8.85  X  10  (farads/meter) 

Relative  permittivity  (dimensionless) 

Intrinsic  impedance  /—  (ohms) 

v  e 

Wavelength  (meters) 

Permeability  a  / n  (henries/meter) 
r  o 

“7 

Permeability  of  free  space  4tt  X  10  (henries/meter) 
Relative  permeability  (dimensionless) 

Reflection  coefficient  (dimensionless) 

Conductivity  (mhos/meter) 

Time  constant  (seconds) 

Azimuth  angle  of  incidence  (dimensionless) 

Elevation  angle  of  incidence  (dimensionless) 

Radian  frequency  2jtf  (second  "S 


56 


Appendix  B 

TRANSMISSION-LINE  FORMULAS 


57 


THIS  PAGE  INTENTIONALLY  BLANK 


TRANSMISSION-LINE  FORMULAS 


Impedance  of  shorted  line  ZQ  tanh  yl  jZQ  tan  $1  ZQ  cog  p &  +  i<yl  sin  pje 


rt, _ +44..  ~  T  .  Approximate  Results  for 

Quantity  General  Line  Ideal  Line 

Low-Loss  Lines 


<D  to  O 
-HOC 
-P  *H  *H 
*H  -P  rH 
-P  *H 
fl  P  bfl 
rt  C  C 
C  etf  O 
cr  3  h 
cf  rt 
TJ 

C  T3  <D 
O  fi  O 
Q>  C 
-P  cd 

3  TJ  -P 

a  a  m 

c  o  -H 
H  J  Q 


O  ^  J  N  ^  I> 
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Wavelength  measured  along  line. 
Velocity  of  light  in  dielectric. 
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CHARACTERISTIC  IMPEDANCE  OF  TRANSMISSION-LINE  CONFIGURATIONS 
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Appendix  C 

CHARACTERISTIC  IMPEDANCE  OF  TRANSMISSION-LINE  CONFIGURATIONS 


T]  e  T1 

=  —  log  cot  —  «  -  log 

7t  S  n 


(0  «  1) 
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CONDUCTIVITY  AND  RELATIVE  PERMEABILITY 
OF  SHIELDING  MATERIALS 
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Appendix  D 


CONDUCTIVITY  AND  RELATIVE  PERMEABILITY 
OF  SHIELDING  MATERIALS 


Material 

Conductivity, 

a 

(mho/m) 

Relative 

Permeability, 

4r 

Silver 

6.8  X  107 

1 

Copper 

7 

5.8  X  10 

1 

Brass 

1.4  X  107 

1 

Aluminum 

7 

3.7  X  10 

1 

Lead 

6 

4.5  X  10 

1 

Nickel 

7 

1.3  X  10 

100 

Mild  steel 

6 

6  X  10 

500 

High-ji  Alloy 

6 

1.6  X  10 

4 

5  X  10 
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Appendix  E 


WIRE  TABLE 


Resistance 

Weight 
(pounds/ 
1000  ft) 

20 Pc 

solid  uopper 

Standard  Annealed — 

Hard- 

100%  Conductivity 

Drawn*^” 

AWG 

Diameter 

Area 

(circular 

mils) 

0hms/1000  ft’ 

* 

2  CPC 

B  &  S 

(mils) 

cPc 

20°  C 

50°  C 

0000 

460,0 

211,600 

0.0451 

0.0490 

0.0548 

0.0504 

640.5 

000 

409.6 

167,806 

0.0569 

0.0618 

0.0691 

0.0635 

507.9 

00 

364.8 

133,077 

0.0717 

0.0779 

0.0871 

0.0801 

402.8 

0 

324.9 

105,535 

0.0906 

0.0983 

0.1099 

0.1010 

319.5 

1 

289.3 

83,693 

0.1141 

0.1239 

0.1385 

0.1273 

253.3 

2 

257.6 

66,371 

0.1440 

0.1563 

0.1748 

0.1606 

200.9 

3 

229.4 

52,635 

0.1815 

0.1970 

0.2203 

0.2025 

159.3 

4 

204.3 

41,741 

0.2289 

0.2485 

0.2778 

0.2554 

126.4 

5 

181.9 

33,102 

0.2886 

0.3133 

0.3503 

0.3220 

100.2 

6 

162.0 

26,251 

0.3640 

0.3951 

0.4418 

0.4061 

79.46 

7 

144.3 

20,818 

0.4590 

0.4982 

0.5570 

0.5120 

63.02 

8 

128.5 

16,510 

0.5787 

0.6282 

0.7024 

0.6456 

49.98 

9 

114.4 

13,093 

0.7297 

0.7921 

0.8856 

0.8141 

39.63 

10 

101.9 

10,383 

0.9203 

0.9989 

1.117 

1.027 

31.43 

11 

90.74 

8,234 

1.161 

1.260 

1.409 

1.295 

24.92 

12 

80.81 

6,530 

1.463 

1.588 

1.775 

1.632 

19.77 

13 

71.96 

5,178 

1.845 

2.003 

2.240 

2.059 

15.68 

14 

64.08 

4,107 

2.326 

2.525 

2.823 

2.595 

12.43 

15 

57.07 

3,257 

2.933 

3.184 

3.560 

3.272 

9.858 

16 

50.82 

2,583 

3.700 

4.016 

4.490 

4.127 

7.818 

17 

45.26 

2,048 

4.665 

5.064 

5.662 

5.204 

6.200 

18 

40.30 

1,624 

5.882 

6.385 

7.139 

6.562 

4,917 

19 

35.89 

1,288 

7.418 

8.051 

9.002 

8.274 

3.899 

20  i 

31.96 

1,022 

9.351 

10.15 

11.35 

10.43 

3.092 
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Appendix  E  (Concluded) 


Solid  Copper 

Resistance 

Weight 
(pounds/ 
1000  ft) 
20°  C 

Standard  Annealed — 

100%  Conductivity 

Hard- 

Drawn*^ 

AWG 

or 

B  &  S 

Diameter 

(mils) 

Area 

(circular 

mils) 

0hms/1000  ft* 

20°  C 

0°C 

20°  C 

5CPC 

21 

28.46 

810.1 

11.79 

12.80 

14.31 

13.16 

2.452 

22 

25.35 

642.4 

14.87 

16.14 

18.05 

16.59 

1.945 

23 

22.57 

509.5 

18.76 

20.36 

22.76 

20.92 

1.542 

24 

20.10 

404.0 

23.65 

25.67 

28.70 

26.38 

1.223 

25 

17.90 

320.4 

29.82 

32.37 

36.19 

33.27 

0.9699 

26 

15.94 

254.1 

37.60 

40.81 

45.63 

41.94 

0.7692 

27 

14.20 

201.5 

47.42 

51.47 

57.55 

52.90 

0.6100 

28 

12.64 

159.8 

59.80 

64.90 

72.57 

66.70 

0.4837 

29 

11.26 

126.7 

75.39 

81.83 

91.49 

84.10 

0.3836 

30 

10.03 

100.5 

95.07 

103.2 

115.4 

106.1 

0.3042 

31 

8.928 

79.70 

119.9 

130.1 

145.5 

133.7 

0.2413 

32 

7.950 

63.21 

151.2 

164.1 

183.5 

168.6 

0.1913 

33 

7.080 

50.13 

190.6 

206.9 

231.3 

212.6 

0.1517 

34 

6.305 

39.75 

240.4 

260.9 

291.7 

268.1 

0.1203 

35 

5.615 

31.52 

303.1 

329.0 

367.9 

338.1 

0.09542 

36 

5.000 

25.00 

382.1 

414.8 

463.8 

426.3 

0.07568 

37 

4.453 

19.83 

481.9 

523.1 

584.9 

537.6 

0.06001 

38 

3.965 

15.72 

607.7 

659.6 

737.5 

677.9 

0.04759 

39 

3.531 

12.47 

766.3 

831.8 

930.0 

854.9 

0.03774 

40 

3.145 

9.89 

966.4 

1049 

1173 

1078 

0.02993 

*  Resistance  at  the  stated  temperature  of  a  wire  whose  length  is  1000  ft  at  20° C. 
^  97.3%  conductivity. 
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Appendix  F 


GROUND  CONDUCTIVITY  IN  THE  UNITED  STATES* 


*  Numbers  on  the  legend,  when  multiplied  by  10  ,  indicate  ground  conductivity 

in  mhos/meter.  (Map  by  FCC.) 
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